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Abstract. The laws of gravitation iiave been tested for a long time with steadily improving precision, 
leading at some moment of time to paradigmatic evolutions. Pursuing this continual effort is of great 
importance for science. In this communication, we focus on Solar System tests of gravity and more precisely 
on possible tests that can be performed with radio science observations (Range and Doppler). After briefly 
reviewing the current tests of gravitation at Solar System scales, we give motivations to continue such 
experiments. 

In order to obtain signature and estimate the amplitude of anomalous signals that could show up in 
radio science observables because of modified gravitational laws, we developed a new software that simulates 
Range/Doppler signals. We present this new tool that simulates radio science observables directly from the 
space-time metric. We apply this tool to the Cassini mission during its cruise from Jupiter to Saturn and de- 
rive constraints on the parameters entering alternative theories of gravity beyond the standard Parametrized 
Post Newtonian theory. 

Keywords: tests of general relativity, radio science 
1 Introduction 

Testing General Relativity (GR) is a long standing and worthy effort in the scientific community. From a 
theoretical point of view, different attempts to quantize gravity or to unify it with other fundamental interactions 
predict deviations from GR. From an observational point of view, cosmological data can not be explained by the 
combination of GR and the standard model of particles, requiring the introduction of Dark Matter and Dark 
Energy. Since these two dark components have not been observed directly, cosmological observations can be a 
hint that the gravitation theory differs from the Einstein theory of gravity. 

Within the solar system, GR is very well confirmed by different types of experiments: tests of the weak 
equivalence principle, Post-Newtonian constraints or fifth-force searches. In Sec. [2l we briefly recall the solar 
system constraints on the gravitation theory and we give motivations to go beyond these stringent constraints. 
In this communication we focus on the possibility to test gravity with radio science measurements (Range and 
Doppler). In order to study the impact of alternative theories of gravity on radio science observables, we present 
a software that simulates Range and Doppler signals directly from the space time metric (and from the initial 
conditions of the bodies considered). This approach allows one to simulate signals in any metric theory. This 
software is presented in Sec. [3] as well as some simulations about the Cassini mission. 
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2 Solar system constraints on gravity 



2.1 Basis of General Relativity 

General Relativity is built on two main principles. The first one, the Equivalence Principle, gives to gravita- 
tion a geometric nature. This principle implies that gravity can be identified to space-time geometry which 
is described by a metric tensor 5^^. Since freely falling test masses follow the geodesies of this metric, their 
motion is independent of their composition. This universality of free fall can be parametrized by a parame- 
ter 77, defined as the relative difference between the accelerations of two test bodies. T he universality of free 
fall has then been tested to an impressive l evel of r? < 10~^^ by Lu nar Laser Ranging (| Williams et al. I I2OO9I) 
and by torsion pendulum (for a review see lAdelberger et al. I I2OO9I) . Let us note that even if this principle 
is one of the best tested in physics some theoretical models coming from unification theories may still pro- 
duce deviations below the current constrai nts (.Damour fc Polvakov 19 941. justifying the necessity to perform 
more accurate experiments like Microscope ( Touboul fc Rodriguesll200ll ) or STEP missions ( Mester et al.l2001 ). 



The Equivalence Principle also defines the behavior of ideal clocks which is independent of their constitution 
and which measures a geometrical quantity t — J y^g^^dxt^ /ds dx'^/ds ds, a point which is also very well tested. 



While the first building block of GR postulates the existence of a metric tensor that determines the trajectory 
of freely falling bodies (the geodesies) and the behavior of clocks, the second building block concerns the form 
of the metric tensor. The metric tensor is determined through field equations, which in the case of GR, are the 
Einstein equations 

_ 1 _ 8ttG 

G^y R^v ^ '^9tii'^ ^4 Tf^iy , (2.1) 

G^n is Einstein curvature tensor, R^^ and R are Ricci and scalar curvatures (these are derived from the metric), 
G is Newton constant, c the speed of light and T^i, the stress-energy tensor. These equations allow one to 
determine the metric tensor from the energy and matter contents of space-time. 

If one considers the solar system and only the Sun's gravitational contribution, the metric tensor can be 
written as an expansion of the Newton potential (pN — (where M is the mass of the central body and r 

a radial coordinate) 

ds^ = {1 + 2(j)N + 2(j)% + . . .)c^dt^ - (1 - 20Ar + . . .)dx^. (2.2) 

This expansion is justified since the value of the Newton potential in the solar system is always smaller than 
(t>N < 10^^. The metric (|2.2p is written in isotropic coordinates. 

The previous two principles produce important effects that need to be taken into account for high precision 
astrometry or for space missions: first, the dynamic is different from Newton gravity (one major consequence 
is the advance of planet perihelia) ; then, the propagation of light is affected ; finally clocks measure a proper 
time which is different from coordinate time t appearing in the metric. All these effects need to be modeled for 
space missions in addition to other classical effects. As a consequence, space missions constitute an excellent 
laboratory to test gravity. 



2.2 Solar system tests of General Relativity 



Gravity tests allow one to compare measurements with the predictions made by some theoretical framework. 
Within the solar system, tw o important frameworks have been developed: the Para metrized Post Newtonia. n 
formalism fully d escribed in Will (1993) and the fifth force search fully described in iTal madge et ah I (Il988l) : 



Adelberger et al. (|2009,) . Both formalisms embed GR metric within a wider class of specifically parametrized 



metric. 



The PPN formalism extends the metric (|2.2p by introducing Post Newtonian Parameters in expansion (|2.2p . 
In the simplest case, two parameters are introduced and the metric describing a spherical and static sun becomes 

ds^ = (1 + 2(j)N + 2/3(t)% + . . .)c^dt^ - (1 - 2-f(j)N + ■ . ■)dx^. (2.3) 

Testing GR with PN formalism consists in measuring the PPN parameters 7 and /? and comparing the exper- 
imental results with their GR value (7 = /? = 1). Thirty years of precise experiments have constrained PPN 



Test of gravity in the Solar Sytem 



239 



parameters very closely around GR (for a review, see IWilll 120061 ). In particular, the observation of the Shapiro 
delay of the Ca ssini probe during a solar conjunction in 2002 gives the best constraint on the 7 PPN parameter 
known to date (jBertotti et al. I I2OO3I) 

7- 1 (2.1 ±2.3) X 10"^ (2.4) 

This constraint is con firmed by deflection measurement with VLBI ( Lambert &: Le Poncin-Lafitte 20091 ) or 
tracking Mars orbiter ([Konopliv et al.ll201l[ ). The present best constraint on the [3 parameter comes from the 
solar system ephemerides (assuming the 7 parameter to be given by (|2.4p ) (jFienga et al.ll201l[ ) 



^- 1 = (-4.1±7.8) X 10"^ (2.5) 
Other type of experi ments also confirmed t his constraint like Lunar Laser Ranging ([Williams et al. I I2OO9I) or 



Mars orbiter tracking ([Konopliv et al.ll201ll) 



The other type of formalism often used in solar system tests of gravity, the fifth-force search, consists in a 
search of the dependence of the Newton potential (j)N with the radial coordinate. This formalism parametrizes 
deviations from the Newtonian potential with a Yukawa potential, justified by unification models. More precisely 
the scale dependance of the Newton potential modifies the temporal part of the metric 



(2.6) 



The a parameter characterizes the amplitude of the deviation with respect to the Newtonian potential and the 
A parameter is a range related to the mass of the new particle that would mediate this fi fth interaction. These 
parameters have been tested in a very wide range (see Fig. 31 of iKonopliv et aL 2011 ). In particular, the a 
parameter is constrained to a very high level of accuracy (a < 10~^°) at Earth-Moon and Sun-Mars distances. 
From this picture, we can also see that windows remain open at very short distances and at very large distances. 



2.3 Is it necessary to go beyond these tests ? 



From the arguments presented in last section, one may wonder if it is necessary to continue to go beyond the 
present constraints on the gravitation theory. The answer is positive for several reasons. Firstly, there exist 
theoretical models predicting deviations smaller than the current constraints. For ex ample, a tensor-scalar 
theory of gravity can naturally be attracted towards GR by a cosmological mechanism (I Pamour Sz Nordtvedt 
I1993I) and lead to a dev iation of the 7 pararneter s maller than the Cassini constraint. Another example is given 
by the chameleon field ( Khourv fc WeltmanI 2004 ) where deviations of GR are hidden in region of high density 
(in the solar system) and remain smaller than current constraints. 

Then, it appears interesting to search for deviations in regions where no test has been performed so far. For 
example, a fifth force may be searched at very small range or at very large range where stringent constraints 
are missing. Let us mention that not all alternative theories of gravity are entering the PPN or fifth-force 
framework. It would be instructive to extend these frameworks to include new types of deviations. To illustrate 
this, we consider two alternative theories of gravity that do not enter the PPN or fifth- force formalism. 

Such a first alternative theory of gravity is provided by Post-Einsteinian Gravity (PEG) ( Revnaud fc Jaekell 
[2005 : Jackcl & Rcynaud 2005, 2006a, b). This theory is based on a non local extension of Einstein field equations, 
as suggested by radiative corrections. Phenomenologically in the solar system, the space time metric can be 
parametrized by two radial dependent potentials S^N{r) and (5$p(r) 



500 

9tj 



= 1 

= 5. 



y (-1 -I- 2</)jv + 25$A, - 25$p) . 



(2.7) 
(2.8) 



This parametrization extends the PPN formalism by replacing parameters by functions. 

Th e second example con sidered in this communication concerns a particular effect due to MOND theory. 
In iBlanchet fc NovakI ( 201ll ) . it is shown that MOND induces within the solar system an External Field Effect. 
This effect is modeled by a quadrupolar contribution to the Newtonian potential 



GM _ 02 
rc? 1(? 



X X-' 



eie, 



(2.9) 



Ci is a unitary vector pointing towards the galactic center and 2.1 10 "^"^ t 
of the quadrupole moment whose value depends on the MOND function. 
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These examples show that there are motivations to improve current tests of gravity but also to look at 
situations previously untested. This can be done by considering theories not entering the traditional frameworks 
or looking in regions of parameters where strong constraints are missing. 



3 Range and Doppler simulations 

In this work we focus on the possibility to perform gravity tests with radio science measurements. This kind of 
test has already been very successful when deriving the PPN 7 constraint with Cassini (see previous section). 
In order to provide estimation of the order of magnitude and of the signature that an alternative theory of grav- 
ity produces on radio science signals (Range and Doppler), we have developed a software that simulates these 
signals. In order to allow for a wide class of alternative theories of gravity, Range and Doppler are simulated 
directly from the space-time metric g^^. This means that it is very easy to change the gravitation theory by 
changing the metric. 



Since GR and the alternative theories of gravity we consider are covariant, we are free to choose the co- 
ordinates to work with (gauge freedom). On the other hand, observations are covariant quantities (or gauge 
independent). Therefore, it is very important to produce simulations that are also covariant. In radio science, 
covariant quantities are based on proper time (time given by ideal clocks). The Range is defined as the differ- 
ence between the reception proper time and the emission proper time and the Doppler is defined as the ratio of 
the received proper frequency over the emitted proper frequency. Three different steps are needed to simulate 
these quantities: the derivation and integration of the equations of motion, the derivation and integration of 
the equations of proper time and the computation of time tra nsfer in curved sp ace-time. These three steps are 
computed from the space-time metric by methods described in Hees et al. ( 201ll) . Finally in order to investigate 



the observables of an alternative theory of gravity in the Range and Doppler data we perform a least-squares 
fit in GR on the different parameters (initial conditions and masses of the bodies) and search for identifiable 
signatures in the residuals. 



In the following, as an example, we present simulations of a two-way Range and Doppler signals for Cassini 
spacecraft from June 2002 during 3 years (when the probe was between Jupiter and Saturn) . A simplified model 
is built with the Sun, the Earth and Cassini spacecraft. We successively consider Post-Einsteinian Gravity and 
MOND External Field Effect. 



3.1 Post-Einsteinian Gravity (PEG) 

As a simplified preliminary study, we focus on the effects of the potential 6^p{r) on the radio science signals. 
Indeed, the potential S(f)pf is already very tightly constrained (from fifth force measurement). We consider a 
series expansion of this potential, that is to say we suppose the spatial part of the metric to be modified as 

m = m\GB. ~ 2% \ Xir + X2r + I (3.1) 

M is the Sun mass, (^7 = 7 — 1 is related to the PPN parameter and xi and xi a-^e parameters characterizing 
deviations due to linear and quadratic terms in the metric. 

We performed different simulations with different values for the three PEG parameters. For example, Fig.[T] 
represents the Range and Doppler residuals due to the presence of a PPN 7 deviation of (S7 = 7 — 1 = 10~^. 
The three peaks occur at solar conjunctions. The signal due to the conjunctions is not a bsorbed at all by the fit 



of the initial conditions. The simulated data are comparable with real data obtained in iBertotti et al.l (|2003l ). 

To summarize the simulations performed. Fig. [2] represents the maximal Doppler residuals due to PEG 
theories of gravity (parametrized by their values of xi, X2 and ^7 = 7 — 1). By requesting the maximal residuals 
to be smaller than the Cassini precision on Doppler signal (lO^^**), we can estimate an upper limit for the PEG 
parameters: xi < 10~^^ , X2 < 2 lO^'^^ wt"^ and ^7 < 3 10^^ (which is very similar to the real estimation 

(Ea). 



3.2 MOND External Field Effect 

In the framework of MOND theory, iBlanchet fc Novak ( 2011 ) have shown that there exists an effect in the 
Solar System, the External Field Effect (EFE), which is modeled by a quadrupole contribution (see Sec. | 
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Fig. 1. Representation of the Range (left) and Doppler (right) signals due to a PEG theory of gravity with 57 = 7 — 1 = 
10~^. The blue (dashed) lines represent the direct difference between simulations in GR and simulations in the alternative 
theory of gravity (with the same parameters). The red (continuous) lines represent the residuals after the least-square 
analysis (the observable signals). 
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Fig. 2. Representation of the maximal Doppler signal due to PEG theory (parametrized by the three parameters xi, X2, 
Sy) for the Cassini mission. The blue (dashed) lines represent the direct difference between the Doppler in alternative 
theories and GR (with the same parameters) while the green (continuous) lines represent the maximal residuals that can 
be observed after the least-square analysis. The red lines represent the Doppler Cassini accuracy. 



Results of our simulations for this effect are given in Fig. [3] where Range and Doppler residuals due to the EFE 
with Q2 = 4.1 10~^^ (the maximal value allowed by the theory) are represented. The residuals are below 
the Cassini accuracy (10~^* in Doppler). Therefore, the Cassini arc considered here appears to be insufficient 
to provide a satisfactory test of MOND EFE. Note that this may be improved by considering other arcs or a 
dedicated space mission. 



4 Conclusions 

Starting from the basis of GR, we have presented current gravity tests performed in the Solar System. We 
have given motivations to increase the accuracy of the current constraints and to look at situations previously 
untested. In this context, we have presented a new software that performs radio science simulations from the 
space-time metric. With this tool, we are able to simulate any space mission in any alternative metric theory 
of gravity and to give the signature of an hypothetical alternative theory of gravity on the Range and/or on 
the Doppler. As examples, we have presented simulations of the Cassini mission in PEG theory of gravity and 
with the MOND External Field Effect. For the PEG theory, constraints have been derived for the parameters 
entering the expression of the spatial part of the metric. For the MOND EFE, we have shown that the predicted 
effect is too small to be detected with the considered arc of the Cassini mission. 
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Fig. 3. Representation of the Range (left) and Doppler (right) signals due to the MOND External Field Effect with 
Q2 = 4.1 10~^^ . The blue (dashed) lines represent the direct difference between simulations in GR and simulations 
in the alternative theory of gravity (with the same parameters). The red (continuous) lines represent the residuals after 
the least-square analysis (the observable signals). 

In the future, further simulations can be done for other theories and other (future and past) space missions to 
answer the following question: can a particular theory of gravity be observed with Range/Doppler measurements 
of a specific mission ? 

A. Hees is research fellow from FRS-FNRS (Belgian Fund for Scientific Research) and he tanks FRS-FNRS for financial support. 
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